Nonlinear conductivity of two-dimensional Coulomb glasses 
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We have studied the nonhnear conductivity of two-dimensional Coulomb glasses. We have used 
a Monte Carlo algorithm to simulate the dynamic of the system under an applied electric field E. 
We found that in the nonlinear regime the site occupancy in the Coulomb gap follows a Fermi- 
Dirac distribution with an effective temperature T^g, higher than the phonon bath temperature T. 
The value of the effective temperature is compatible with that obtained for slow modes from the 
generalized fluctuation-dissipation theorem. The nonlinear conductivity for a given electric field and 
T is fairly similar to the linear conductivity at the corresponding Tcff. We found that the dissipated 
power and the effective temperature are related by an expression of the form {T^g — T°')T^g°' . 

PACS numbers: PACS numbers: 72.20.Ht, 72.20.Ee, 72.80.Ng 



I. INTRODUCTION 

Electron transport in Coulomb glasses (CG) has been 
investigated for decades. These glasses are systems with 
electronic states localized by the disorder and long range 
Coulomb interactions between carriers. At low tempera- 
tures, conductivity in CG is by hopping, where the transi- 
tion rates for electron jumps depend exponentially on an 
energy factor and on a spatial factor. In the limit of very 
low temperatures, an optimization of the total penalty 
paid through these two factors leads to the mechanism 
coined by Mott as variable range hopping (VRH). Mott 
obtained the precise law for the DC conductivity a of 
non-interacting systems in this regimei and Efros and 
Shklovskii (ES)^ modified his argument to include the 
effects of Coulomb interactions by considering the spe- 
cific form of the single-particle density of states in CG. 
The conductivity in this case is of the form 



a cx exp{- (To/T)'/"} 



(1) 



with the exponent 1/2 independent of the dimensionality 
of the system. Tq — jSe^ /{ek^) is a characteristic temper- 
ature, e the dielectric constant of the material, k Boltz- 
mann constant, ^ the localization radius of the electrons, 
and /? a numerical coefficient that depends on dimen- 
sionality. From percolation theory for hopping transport 
one obtains /3 ~ 6.5 for 2D systems (see references inQ). 
The conductivity of many different types of systems have 
been found to obey this law, Eq. ([1]). 

Nonlinear effects in electron transport are specially im- 
portant in CG, where interactions usually increase non- 
linearities and can also establish an effective temperature 
for the electrons, higher than the phonon bath temper- 
ature, at very low temperatures when the thermal con- 
ductance between the electrons and the phonons is not 
large enough to dissipate all the electrical power pro- 
vided to the system. While there are many experimental 
studies of nonlinear effects on systems showing VRH in 
the presence of a Coulomb gap, there is no proper the- 
ory addressing the problem. Some experimental results 
have been interpreted in terms of theoretical approaches 
designed for the noninter acting VRH regime, and ex- 



tended to the interacting case by adapting the typical 
hopping length. For electric fields satisfying the condi- 
tion Ec < E < kT/elo, where Ec — kT/eL, Iq is the aver- 
age distance between nearest neighbor impurities and L 
is some typical hopping length, the theoretical modela^*^ 
predict a conductivity of the form 



cr(T, E) = (t{T, 0) exp {eEL/kT} . 



(2) 



The different models differ in the precise form of the typ- 
ical hopping length L. While the argument of PoUak and 
RiessA adapted to VRH in CG would produce L cx T"^/^, 
that of Shklovskii^ would result in L oc T~^. These ap- 
proaches were called field-effect models and where rela- 
tively successful in explaining early experimental results. 
Some recent experiments can also be interpreted in terms 
of these models, for example, the results of Grannan et 
al.^, which fit well to L oc T^^, and those of Zhang et 
al^, which fit adequately to L c>c T^^/^. 

In 1990, Wang et aZ*§. presented thorough results on 
neutron-transmutation doped Ge that could not be fit- 
ted with any field-effect model and which could be ex- 
plained with the so called hot-electron model. This model 
assumes that the applied electric power is deposited in 
the electron system and is dissipated to the phonon sys- 
tem through the electron-phonon coupling. For slowly 
relaxing systems in general, and hopping systems at very 
low temperatures in particular, this coupling is too weak 
to be able to thermalize together the electron and the 
phonon systems. The electron-electron interaction estab- 
lishes an effective temperature T^g for the electron sys- 
tem, which remains higher than the phonon bath temper- 
ature T. The model also assumes that the conductivity 
for any field E and temperature T only depends on the 
electron system and is the same as the linear conductivity 
at the effective temperature 

a(r, E) = a(reff, 0) cx exp{- {To/T.sf^}. (3) 

The effective temperature can be calculated from the 
thermal conductivity between the conducting electrons 
and the phonon bath. This quantity is not known for 
hopping systems and by analogy with metals one assumes 
that it is proportional to a power of the temperature. 
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Then the effective temperature is related to the electric 
power P through 

P^a{Tl^^-T^) (4) 

where a is independent of temperature. 

Zhang et ali- performed a systematic study of non lin- 
ear effects on doped Si and Ge and concluded that field- 
effect models explain adequately the experimental results 
when Tq/T > 135, while hot-electron models provide a 
better fit when Tq/T < 135. In line with these results, 
Gershenson et al^ studied two-dimensional hopping and 
concluded that in systems with small localization lengths 
(and so large Tq) field-effects dominate, while in systems 
with large localization lengths the hot-electron model ex- 
plains adequately the results. Recently, many exper- 
imental results have been fitted with the hot-electron 
modeli2r2ij while some still have been interpreted in 
terms of field-effect models^^'^^. 

The hot-electron model together with a strong temper- 
ature dependence of the conductivity presents an insta- 
bility at very low temperatures. As the applied voltage 
is increased the current may abruptly change by several 
orders of magnitude, as it has been observed in amor- 
phous indium oxide films in the insulating state^i. This 
effect has been explained in terms of a bistability of the 
effective temperature^^. We will point out later on that 
a similar effect could also be present in Coulomb glasses 
at sufficiently low temperatureai^. 

In this paper, we perform Monte Carlo (MC) simula- 
tions of CG in the nonlinear regime and calculate the 
conductivity as a function of temperature and applied 
field. As we can measure the effective temperature of the 
electrons directly^ and the absorbed and emitted powers 
separately, our results are an excellent tool to check the 
validity of the hot-electron model. After presenting our 
model in the next section, we study the effective tempera- 
ture Toff as obtained from the site occupation and analyze 
its relation with the temperature of slow modes obtained 
from the generalization of the fluctuation-dissipation the- 
orem. In section IV, we present the results for the nonlin- 
ear conductivity in the variable range hopping regime as 
a function of the temperature and the electric field and 
analyze the applicability of the hot-electron model. In 
section V, we study the dissipated power in the non lin- 
ear regime in terms of Teg. In section VI, we analyze the 
predictions of the hot-electron model. We finalize with 
some conclusions. 

II. MODEL 

To calculate the conductivity of two-dimensional 
Coulomb glasses we consider the HamiltonianSi 

(0, -f x,E)n, + — — . (5) 

where = 0, 1 are occupation numbers, K is the com- 
pensation, equal to 1/2, and are the random site 



energies chosen from a box distribution with interval 
[— VF/2, W^/2]. E is the value of the applied electric field. 
Xi is the coordinate along the direction of the applied 
field of site i and is the distance between sites i and 
J . We consider square samples of lateral size L and with 
N sites placed at random with a minimum separation 
between them of 0.2Zo- = L/^/N is our unit of dis- 
tance and 1/Zo is our unit of energy and temperature. 
We study systems ranging from 500 to 4000 sites with a 
range of disorder W = 2 and two values of the localiza- 
tion length, ^ = 1 and 2. We implement periodic bound- 
ary conditions and Coulomb interactions are calculated 
using the minimum image convention^^i, which is a rea- 
sonable choice for disordered systems in order to reduce 
finite size effects avoiding at the same time artificial long- 
range correlations. The electric field energy increases lin- 
early in the direction of the applied field. Once periodic 
boundary conditions are imposed, the total energy of the 
system is not well-defined, but the transition energies are 
well-defined provided that we restrict ourselves to hops 
of length shorter than half the lateral size of the system. 
The previous arrangement induces a permanent current 
in the direction of the field when we set up our dynamical 
procedure. 

We concentrate here in the regime where the tempera- 
ture is low enough for conduction to be by variable range 
hopping, and high enough so that a stationary state is 
achieved relatively fast. In this case we can consider only 
single-electron transitions, since we are not in the regime 
where many-electron transitions are importanlj^^. The 
clusterization algorithm previously developed by us^^ is 
not efficient in the temperature range considered here. As 
we are studying stationary states, not equilibrium prop- 
erties, a kinetic MC algorithm that keep track of the 
physical time taken in the different processes must be 
used^^. The MC method employed first chooses a pair of 
sites with a probability proportional to exp(— 2rij/^}^. 
Then if one of the sites chosen is occupied and the other 
empty it exchanges their occupation when the total tran- 
sition energy AE, including the energy due to the applied 
field, is negative or with probability exp(— AE'/T) when 
AE > 0. The time step of our MC procedure is then 
equal to Tq/ exp(— 2rij/^), where tq is the inverse 
phonon frequency, of the order of 10"^"^ s2L2^. 

In the non linear regime we define the conductance as 
the ratio of the intensity divided by the applied voltage, 
V = EL. As we consider square samples, the conduc- 
tivity is thus a = I /{EL). We start from a random 
configuration and follow the dynamics at a given tem- 
perature, monitoring all relevant magnitudes. Once we 
are in a stationary situation, we obtain the conductivity 
of each sample through the displacement of the center 
of mass of the electrons for a given time interval^S. The 
temperature range studied goes from 0.05 to 0.2, the typ- 
ical number of samples is 5000 and the simulation time 
is IO^tq. The values of the applied electric field run up 
to 2T. 
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III. EFFECTIVE TEMPERATURE 

With our MC procedure we simulate the dynamics of 
CG in the presence of a relatively strong electric field, of 
the order of T, so that we are in the nonlinear regime. 
Once the system has reached a stationary state, we cal- 
culate the site occupation probability /(e^) as a function 
of site energy Cj, which is defined as 



(rij - K) 



(6) 



We note that this energy does not include the electric 
field contribution. We find that /(ei) follows pretty well 
a Fermi-Dirac distribution with an effective temperature 
Teff, which depends on both T and the applied field E. 
In Fig. f we represent the occupation probability of sites 
near the chemical potential as a function of site energy for 
T = O.f and two different values of the electric field E = 
r/2 (circles) and T (squares). The continuous lines are 
Fermi-Dirac distributions corresponding to Teg — 0.112 
and 0.136, respectively. 




FIG. 1: (Color online) Occupation distribution function near 
the chemical potential at T = 0.1 for two values of the ap- 
plied field E = r/2 (blue circles) and E = T (red squares). 
The continuous lines are Fermi-Dirac distributions for tem- 
peratures 0.112 and 0.136, respectively. 

Let us call n{AE) the density of electron-hole exci- 
tations of energy AE and longer than a given length, 
that we take equal to 10. The relative probability of 
having long electron-hole excitations with negative and 
positive energies, n{— AE) / {n{— AE) + n{AE)), follows 
a distribution very similar to the site energy occupation 
probability. The situation is similar to simulations of 
relaxation^*^. Long jumps and injection or extraction of 
single particles near the Fermi Level are both very slow 
process equilibrated at the same effective temperature. 
In practice, it is easier to calculate Tos from the distri- 
bution of the excitations than from the site occupation 
function, since for the former the statistics are better and 



the procedure is insensitive to the determination of the 
chemical potentiaP^. 

The quality of the fit of the data in Fig. 1 to a FD 
distribution indicates that the electron-electron interac- 
tion has indeed thermalized the electrons near the chem- 
ical potential at the effective temperature. However, we 
note that not all electrons are thermalized at T^s- Fast 
electronic modes, like short electron-hole excitations, are 
equilibrated at T, while slow modes, like long excitations, 
are equilibrated at Toff. 

The existence of two different temperatures in the 
electronic system can be also observed through the ex- 
tension of the fluctuation-dissipation theorem for non- 
equilibrium systems in the way described in the spin glass 
literature^i^. In order to compare our results with those 
based on the fluctuation-dissipation theorem, we follow 
the same procedure as in Refs. |32||33| . We study the re- 
sponse of the system to an external perturbation of the 
form Stpi — \{t)ipi = XoLpi0{t — tw), where Aq ^ 1 and ipi 
are normalized random variables uncorrelated from site 
to site and from the original random site energies. The 
quantity conjugated to X{t) is dn{t) = 1/N ^^{Sni{t)(fi) , 
where Sni{t) = ni{t) — K. In the linear response regime 
the latter quantity is given by 



Sn{t) = Ax(t -I- t^,t^), 



(7) 



where x is the local susceptibility. At equilibrium, x is 
related to the local charge correlation function 



C(t -t- tyj, t^) 



N 



^{5n,{t + t^)5n^{ty,)), (8) 



through the fiuctuation-dissipation theorem Tx{t) = 1 — 
C (t) . Out of equilibrium, this relation between x and C 
is not verified in general, but in slowly relaxing systems 
it is often found that at short time scales the equilibrium 
relation still holds, while at long time scales there is also 
a linear relation with a different slope. From this slope 
we can define an effective temperature through dx/dC — 

-1/Tefi^. 

We calculate C{t + tw,tw) and x{t + tw,tw) for several 
values of the applied electric field and in Fig. [5] we present 
parametric plots of Tx{t -I- in, , in, ) versus C{t + t.u],tw) for 
T — 0.1. The system is out of equilibrium in a stationary 
state due to an applied electric field, and so t^j is irrel- 
evant in this case. Each curve in Fig. [2] corresponds to 
a different electric field E — 0.008, 0.1 and 0.15, from 
top to bottom. All curves start, for short time intervals, 
on the lower right corner with an slope equal one and 
then fiatten out on the left part, for long time intervals, 
presenting a slope smaller than one and so an effective 
temperature higher than T. The dashed lines have slopes 
equal to T/Tcff, where the effective temperature has been 
calculated from the excitation probabilities as described 
above. 

In Fig. [5] we can appreciate a fairly well defined effec- 
tive temperature for slow processes, much better defined 
than in the case of relaxation from high energy states^. 
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FIG. 2: (Color online) Tx{t + tw,tw) as a function of the 
correlation function for E = 0.008 0.1 and 0.15, from top to 
bottom, and T — 0.1. The dashed lines have slopes equal to 
—T/Teft as obtained from the excitation probabilities. 



We also note the good agreement between the Tcs ob- 
tained from the long excitation probabilities and from 
the extension of the fluctuation-dissipation theorem. We 
would like to add that the use of the latter approach is 
a noisy procedure that requires averages over many sam- 
ples, much less convenient than our method, although 
this needs the explicit identification of slow processes. In 
spin glasses, for example, we have not been able to apply 
our method so far. Coulomb glasses are very suitable to 
study the role of an effective temperature since we have 
two types of simple slow processes: the changes in oc- 
cupation of sites near the chemical potential and long 
electron-hole excitations^S.. 

We have systematically studied the effective tempera- 
ture in terms of the electric field and temperature. From 
now on we will refer to the effective temperature ob- 
tained through the density of excitations, which is the 
most precise and less noisy procedure. In Fig. [3] we plot 
Toff as a function of E for several values of the temper- 
ature T = 0.06 (circles), 0.08 (up triangles), 0.1 (down 
triangles) and 0.12 (diamonds). We have also considered 
two values of the localization length ^. Solid symbols cor- 
respond to ^ = 1 and empty symbols to ^ = 2. The lines 
are a guide to the eye. These data cannot be obtained 
experimentally, thus our direct calculation of Tcff is a 
powerful tool to evaluate the validity of the hot-electron 
model. After we study the non linear conductance, we 
will come back to this issue. 

In the presence of an electric field, the electronic sys- 
tem can increase its energy by long hops along the direc- 
tion of the field. Part of this energy is very quickly dis- 
sipated via phonon emission by fast degrees of freedom, 
like soft dipoles and possible relaxation of the electron 
that performed the long jump. Due to the Coulomb in- 
teraction, the energy provided by the field is distributed 
over the whole electronic system and part of it excites 
slow degrees of freedom. The system (more specifically. 




FIG. 3: (Color online) Effective temperature as a function of 
the applied electric field for several temperatures T = 0.06 
(red circles), 0.08 (green up triangles), 0.1 (blue down trian- 
gles) and 0.12 (magenta diamonds). Solid and empty symbols 
correspond to ^ = 1 and 2, respectively. 

its slow degrees of freedom) cannot relax fast enough and 
remains hot. In the absence of slow degrees of freedom, 
one would expect an effective temperature TcS = 
(for large enough values of the electric field) as derived 
in Ref. The linear dependence between Tes and E 
at high fields in Fig. [3] may correspond to this situation. 
The slopes of these straight segments for ^ = 1 are indeed 
very close to 0.5, while the slopes for ^ = 2 are 0.95 in 
agreement with this theoretical prediction. This model 
cannot explain the behavior for small fields, in particu- 
lar the region for intermediate fields where Tcff is almost 
independent on ^, and the change in curvature for low 
fields and low temperatures. 

In the nonlinear regime, the electric field produces a 
partial filling of the Coulomb gap. This effect can be 
explained in terms of the effective temperature. In equi- 
librium, the density of states at the Fermi Level, g(0), 
depends on T and in 2D this dependence is linear at 
low temperatures. We previously showed that in relax- 
ation experiments g{Q) as a function of Tcff presents the 
same dependence as in equilibrium^*'. Here we found that 
the same relationship holds in the non linear regime. In 
conclusion, the quantity g{0), which could be measured 
experimentally, constitutes an excellent thermometer for 
the electrons effective temperature. 



IV. NONLINEAR CONDUCTIVITY 

We have calculated the conductivity for different values 
of the applied field and the temperature in the variable 
range hopping regime. We observe a very small linear 
regime for small values of the field and then a systematic 
increase of the conductivity for high values of E. This 
increase also depends on T . In Fig. 3] we represent the 
conductivity cr(T, E) on a logarithmic scale as a function 
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of E for the same values of T as in Fig. [3l The tem- 
peratures are T — 0.06 (circles), 0.08 (up triangles), 0.1 
(down triangles) and 0.12 (diamonds), and the localiza- 
tions lengths are ^ = 1 (solid symbols) and 2 (empty 
symbols). The lines are a guide to the eye. For high field 
values, the temperature dependence of a{T,E) is weak. 
For our model, the conductivity tends to zero for very 
high values of E^^ (not shown in the figure) . Later on we 
will analyze how to collapse on a single curve the data for 
different values of T and the applicability of the so-called 
field-effect models and the hot-clcctron model. 




FIG. 4: (Color online) Conductivity on a logarithmic scale 
as a function of E for T — 0.06 (red circles), 0.08 (green 
up triangles), 0.1 (blue down triangles) and 0.12 (magenta 
diamonds), and ^ = 1 (solid symbols) and 2 (empty symbols). 

The conductivity in the linear regime follows the T~^/^ 
law as can be seen in Fig. [SI where cr(T, 0) is represented 
by large squares on a logarithmic scale as a function of 
rp-i/2^ The lower set of data (solid symbols) correspond 
to ^ = 1 and the upper set (empty symbols) to ^ = 2. 
The straight lines are fit to these linear regime data. The 
relatively good fitting of the data to a straight line indi- 
cates that the standard form for variable range hopping 
in the presence of interactions, Eq. (IT|), is fairly well sat- 
isfied. An slightly better linear fit is obtained if we rep- 
resent cr(T, 0)T, instead of cr(T, 0), avoiding certain log- 
arithmic corrections. The characteristic temperature Tq 
would also be in that case in better agreement with previ- 
ous results^. However, we prefer to represent a{T, 0) it- 
self to keep the strongest possible analogy with the usual 
analysis of experimental results in this area and to avoid 
having to rescale the vertical axis with Teff. Our conclu- 
sions do not depend on this choice. From the slope of the 
straight lines we found that the characteristic tempera- 
ture To is in this case equal to 4.43 for ^ = 1 and 2.62 
for ^ = 2. We found that the constant of proportionality 
implicit in Eq. ([T]) slightly depends on ^. 

As we can determine Tcfi directly from the long one- 
electron excitation probabilities, we can check the va- 
lidity of the hot-electron model, which assumes that 
at a given E and T the conductivity is the same as 




FIG. 5: (Color online) a-{T, 0) on a logarithmic scale as a 
function of T^^^'^ (large squares) for ^ = 1 (solid symbols, 
lower set of points) and 2 (empty symbols, upper set). We 
also plot (T(r, E) in the nonlinear regime as a function this 
time of T^g . The data and the symbols employed in the 
non linear regime are the same as in Fig. 3] 

the linear conductivity at the corresponding Teg, i.e., 
(j{T,E) = tT(Teff,0). In Fig. [5J we have also represented 
cr(T, E) on a logarithmic scale as a function of the corre- 

— 1/2 

spending T^g for T = 0.06 (circles), 0.08 (up triangles), 
0.1 (down triangles) and 0.12 (diamonds). The solid sym- 
bols correspond to ^ = 1 and the empty symbols to 2. 
We see that the data in the nonlinear regime lie rela- 
tively close to the T"^/^ law corresponding to the linear 
regime. The difference between a{T, E) and a{TcS: 0) 
is at most a factor of two. The quality of the agree- 
ment of the data for ^ — 2 with the hot-electron model 
is quite good, much better than for the case ^ = 1. This 
trend is in line with experimental results that show that 
the hot-electron model works better for lower values of 
To/T, which in our simulations is equivalent to larger val- 
ues of ^. The increase of cr as a function of the electric 
field can be roughly interpreted as due to an increase in 
Toff. The energy factor entering in the transition rates 
between two configurations can be decomposed into two 
terms, one due to the occupation of the configurations 
and the other to the energy penalty of the transition. 
The first one depends on the effective temperature, while 
the second one depends, in principle, on the phonon tem- 
perature T. Fig. [S] indicates that the contribution of the 
first term masks to a large degree the effect of the second 
term. 

In order to further check the predictions of the hot- 
electron model we must study the power dissipated by 
the system. 



V. DISSIPATED POWER 

The relation between the effective temperature and the 
electric field must be linked to the energy dissipation ca- 
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pability of the syste m^^'^^ . Most authors have assumed 
that the power supplied by the electric field is related 
to Tcff through Eq. (|4]). This type of expression was ob- 
tained for metals, and it was also used, with different val- 
ues of the exponent /?, in the the diffusive regime and even 
in the localized regime without any theoretical support. 
Tcff is usually obtained from the hot-electron assumption, 
that is, through the relation a{T,E) = cr(roff,0). In the 
cases where the hot-electron model works well, the power 
provided by the electric field is indeed approximately re- 
lated to the effective temperature by Eq. (|4]). 

In our simulations, we can obtain more detailed infor- 
mation than in experiments since we can measure Tcff 
directly from the excitation spectra, and because we can 
compute separately the powers absorbed and emitted by 
the electrons. 

We have first calculated the power absorbed by the 
electronic system from the phonons at equilibrium at a 
temperature T in the absence of an electric field (in this 
case, this power is of course equal to the power emitted 
back by the electrons). In Fig. [6] we show the absorbed 
power in equilibrium as a function of T. This power has 
been obtained by two different methods: Monte Carlo 
simulations (solid symbols) and through a set of low 
energy configurations (empty symbols). Different sym- 
bols correspond to different system sizes. In the range 
0.001 < T < 0.3, wc found that the data can be fitted 
by bTf, with 7 = 2.151 ± 0.001, independent of the lo- 
calization length ^, and b a constant weakly dependent 
on ^. The value of 7 is much smaller than the values of 
the exponents found in the expressions for the dissipated 
power, usually between 5 and 6 (see, for example, Ref. 
0). This discrepancy is due to the fact that most of the 
exchange of energy in equilibrium is absorbed and emit- 
ted by short excitations, contributing as T^, while most 
of the dissipated power in DC current is related to more 
complex excitations. 

We now apply an electric field and measure the ab- 
sorbed and emitted powers in terms of T and E. The 
difference between them is, of course, the electric power 
provided by the field a{T,E)E'^. We observe that the 
power emitted by the system is equal to the value at 
equilibrium plus a (usually smaller) term that only de- 
pends on Toff and can be pretty well represented by a 
power law, with an exponent close to 5. However, the 
absorbed power is equal to the power in equilibrium plus 
a term that depends on both T and Tcff. This is in dis- 
agreement with the standard assumption, i.e., that it is a 
function of T only. We found that the power provided by 
the electric field can be fitted pretty well by an expression 
of the form 

P = (7(T, E)E'^ = c (T<fff - T") Tfg-". (9) 

The best fit corresponds to a = 1.99 and /3 = 5.26. We 
have considered a = 2 and (3 — b for simplicity and be- 
cause the fit is also good for these exponents. We note 
that the right-hand side in Eq. ([9]) is the difference be- 
tween the powers emitted and absorbed by the electrons. 




T 



FIG. 6: (Color online) Absorbed power in equilibrium as a 
function of T obtained with Monte Carlo simulations (solid 
symbols) and through a set of low energy configurations 
(empty symbols). Different symbols correspond to different 
system sizes: A'^ = 500 (green triangles), 1000 (red circles) 
and 2000 (black squares). 



but each term separately does not correspond to the full 
emitted and absorbed powers, since the large contribu- 
tion by short dipoles is equal for both powers and can- 
cels. In Fig. [3 we represent (T^^g - T^) T^^^ as a func- 
tion of the dissipated power P for several temperatures, 
T — 0.06 (circles), 0.08 (up triangles), 0.1 (down trian- 
gles) and 0.12 (diamonds). Solid symbols correspond to 
a localization length ^ = 1 and empty symbols to ^ = 2. 
The straight lines are linear fit to the data, forced to 
pass through the origin, for ^ = 1 (continuous line) and 
2 (dashed line). We note both the good overlap of the 
data and the global linear behavior. The proportionality 
constant c depends strongly on ^. 

An alternative and more demanding way to represent 
the data in figure [3 usually employed in the literature^, 
is by plotting the electric power plus the T dependent 
contribution in Eq. ^ (the extra absorbed power) as a 
function of Teff. In Fig. H we plot a(T,E)E'^ + cT^T^g, 
where c is the inverse of the slope of the straight line in 
Fig. El versus Tgff on a double logarithmic scale. We again 
see a good overlap of the data and a fairly linear behavior, 
indicating in this case a good power law dependence of 
the emitted power on Toff. The straight line is a linear 
fit to the data and its slope is {3 — 5.25. This value of (3 
is the same as our best estimate for /3 obtained from an 
overall nonlinear fitting, which is not surprising since the 
increase in the emitted power with respect to equilibrium 
is much larger than the increase in the absorbed power. 

In the inset of Fig. [5] we represent a(T, E)E'^ + cT^ as 
a function of Toff for the same data as in the main part 
of the figure. It clearly shows that the data for small 
electric fields deviate from a straight line and do not 
overlap for different temperatures. A similar behavior is 
observed in several experimental results*. Thus, we con- 
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FIG. 7: (Color online) [T'^g - T^) T^g versus P for T = 0.06 
(red circles), 0.08 (green up triangles), 0.1 (blue down trian- 
gles) and 0.12 (magenta diamonds), and ^ = 1 (solid sym- 
bols) and 2 (empty symbols). The straight lines are a linear 
fit passing through the origin for ^ = 1 (continuous) and 2 
(dashed). 
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FIG. 8: (Color online) a{T,E)E^ + cT'^T^g as a function of 
Tcff for the same data as in Fig. [T] Inset: (t(T, E)E'^ + cT'' as 
a function of Tofi for the same data as in the main part. 



elude that Eq. ^ for the dissipated power is better than 
Eq. (U), usually employed in the literature. The use of 
Eq. ([9]), instead of Eq. (|4]), can explain many "anomalies" 
found in experiments for low dissipated powers as, for 
example, in two-dimensional interface a^"'^^ , in neutron- 
transmutation-doped GeS. and in amorphous indium ox- 
ide films in the insulating regime^^. 

The values of the exponents a and /3, entering in Eq. 
found in our simulations depend on the model used 
and it may be different in experimental situations, de- 
pending for example on the density of phonons and on 
the dimensionality of the system. The corresponding ex- 
ponents may depend on the particular situation studied 
and remain to be measured experimentally or properly 
calculated. The important message is that the absorbed 



power depends on both T and T^s and the dissipated 
power is given by Eq. © , instead of Eq. ^ . 



VI. HOT-ELECTRON MODEL 

We have combined Eqs. ([3]) and ([9]) to obtained our 
modified version of the hot-electron model. The value of 
Tq and the proportionality constant in Eq. ([3]) are empiri- 
cally obtained from the fit of the data in the linear regime 
in Fig. [SJ In Eq. © we take a = 2 and /3 = 5 and the 
constant c is obtained from Fig. [71 With these values of 
the parameters we find numerically Toff and cr(T', E) for 
each value of T and E. We note that there are many dif- 
ferent implementations of the hot-electron model in the 
literature and the one more often found in experimental 
papers is not the same as ours. The latter obtains Toff 
empirically from the relation a{T,E) — a{Tcff,0) and 
uses this value to fit the expression of the power. 

We found empirically that the values of cr(T, E) for dif- 
ferent temperatures and electric fields can be overlapped 
relatively well by plotting a{T, E)/a{T,0) as a function 
of E/T^. In Fig. 1 we plot a{T,E)/a{T,0) (symbols) 
on a logarithmic scale versus E/T'^ for the same data as 
in Fig. m The curves correspond to the predictions of 
the hot-electron model for ^ = 1 (continuous line) and 
2 (dashed line) and T = 0.1. The predictions for other 
values of T overlap between themselves relatively well for 
each ^. The agreement between the hot-electron model 
and our simulations is reasonable, specially for ^ = 2. 




5 10 15 20 

E/f 



FIG. 9: (Color online) cr(r, E)/a{T, 0), on a logarithmic scale, 
as a function of E/T'^ for the same data as in Fig. 3] The 
curves correspond to the predictions of the hot-electron model 
for 5=1 (continuous line) and 2 (dashed line). 

According to field effect models, Eq. ©, 
\og{a{T,E)/a{T,Q)) should be proportional to EL/T. 
We found that this dependence is not fully verified, but 
there is no strong disagreement either. The collapse 
of the conductivity data as a function of E/T^ would 
indicate that the characteristic hopping length for 
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nonlinear effects would be inversely proportional to T 
as predicted by Shklovskii's model. However we note 
that the localization length dependence is not the one 
expected. 

VII. CONCLUSIONS 

We observed that slow electronic modes of CG are 
thermalized at an effective temperature higher than T 
in the non linear regime. Although this effective tem- 
perature is well defined, as we verified through the occu- 
pation distribution, the density of long excitations and 
the fluctuation-dissipation theorem, the detail mecha- 
nism underlying this thermalization is not clear and de- 
serves farther investigation. For moderate electric fields, 
the conductivity increases with the field and this incre- 
ment can be roughly interpreted in terms of the increase 
in the effective temperature. 

The non linear conductivity in terms of the effective 
temperature roughly follows the T^^/^ law with a char- 
acteristic temperature Tb similar to that of the linear 
regime. From a theoretical analysis one could conclude 
that this To should be different from that of the linear 
regime and even may depend on T. According to our 
results these deviations from the linear regime are small. 

The power provided by the electric field in the strongly 
localized systems that we consider can be expressed as a 
difference of powers of Teg and T . However, the expres- 
sion that fits better the numerical results, Eq. ([9]), is not 



the one usually employed in the literature. It would be 
interesting to find the exponent a and /3 from experimen- 
tal results. 

The switching transition observed in an amorphous in- 
dium oxide film in the insulating state^^ was explained 
in terms of a bistability of the electron temperature Teg, 
as obtained from the hot-electron model^S'. This effect 
might be observed in CG in the variable-range hopping 
regime with Tg /T ^ 1 or in another hopping regime with 
a stronger temperature dependence of the conductivity. 
The possible mechanism for switching in our system is re- 
lated to the slow degrees of freedom. The current heats 
the electronic system and part of this heat remains in the 
slow modes producing an effective temperature, which 
strongly increases conductivity. In fact Eq. ([9|) for the 
power, with a = 2 and /3 = 5, and Eq. ^ for the con- 
ductivity, predict that the bistability will occur at tem- 
peratures lower than To/(144-\/3) and in this case the 
hot-electron model should work pretty well. This tem- 
perature range is difficult to simulate with our present 
MC method due to the glassy nature of the system, but 
this is an interesting problem that deserves further study. 
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